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Preface, 


The  purpose  of  the  study  was  to  develop  an  extremely  user  friendly 
computer  code  to  simulate  an  injection  locked  pulsed  C02  laser.  The  use 
of  this  program  is  most  suited  to  the  design  of  lidar  applications. 

The  program  was  written  under  the  BASIC  language  format  for  use  on 
an  IBM-PC/AT  compatible.  Unfortunately  a  version  of  BASIC,  Microsoft's 
BASIC  version  7.0,  was  not  available  until  after  most  of  the  code  was 
written.  Microsoft's  BASIC  will  allow  program  and  data  modules  to 
reside  in  extended  or  expanded  memory.  This  advantage  would  result  in 
more  longitudinal  modes  tracked,  longer  pulse  times,  and  smaller 
interval  sizes  used  in  computations. 

I  wish  to  thank  MAJ  Stone  for  his  efforts  and  understanding  during 
my  thesis. 

Allen  M.  Susie 
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The  purpose  of  this  chapter  is  to  define  key  terms  and  to  detail 
the  direction  and  scope  of  the  work.  This  chapter  is  divided  into  three 
sections.  The  first  section  consists  of  the  definitions  used  throughout 
the  document.  The  second  section  details  the  purpose  of  the  work 
effohrt.  The  final  section  states  the  scope  of  the  work. 


With  the  discovery  of  the  laser  there  have  toe'dii  continual  attempts 


to  increase  the  "spectral  brightness'XOTl  :229)T  of  a  laser  for 
applications  such  as  signal  processing^-and^lidar .  To  this  end 
techniques  using  either  passive  or  active  cavity  elements  have  been 
developed  to  isolate  a  single  loqgAfcudinal  mode  in  both  continuous  and 


pulsed  laser  systems  (ffl : 230)5  .  Although  passive  elements  have  useful 
applications  in  continuous  systems  their  use  in  pulsed  systems  is 
difficult  to  control  (^1x2  .  The,, An jec t i on  of  an  external  signal  has 

proved  very  successful  /ll : 230%  Injection  jiode  locking  is  a  technique 
where  a  single  frequency  signal  is  usedjto  seed  the  active  medium  during 


the  gain  build-up  period  (11:238).  The  laser  transitions  in  C02  gas 
lasers  are  due  to  changes  in  the  molecular  energy  in  the  vibraticnai  or 


rotational  aspects  of  the  molecule.  C02  may  transition  in  one  cf  three 
bands;  the  Regular  band,  the  Hot  band,  and  the  Sequence  band.  In  the 
Regular  band  transitions  the  energy  exchange  between  vibrational  states 
is  very  rapid  because  of  the  near-resonant  energy  transfer  (13:68).  The 


energy  in  one  aspect  of  excitation  can  be  exchanged  for  another  mode  of 
excitation.  For  the  purposes  of  this  thesis  the  focus  will  be  on  the 
transition  between  the  first  asymmetric  stretch  mode  and  the  first 


1 


symmetric  stretch  mode  in  the  Regular  Band  (13:14)  or  the  00°1  to  10°C 
10.6  urn  transition. 

Due  to  collisions  between  molecules  and  the  effect  of  the  movement 
of  molecules  inside  of  the  cavity  the  stimulated  emission  cross  section 
becomes  broadened.  At  low  pressures  below  10  torr  the  effect  is 
primarily  due  to  Doppler  shifts  of  the  transition  frequencies  (13:59). 
Above  50  torr  the  collision  processes  dominate  (12:147;  13:60).  This  is 
called  collision  broadening  or  pressure  broadening.  This  model  is 
limited  to  the  pressure  broadened  regime. 

P.urpo.se  .of  JtheJforls 

An  injection  locked  laser  delivers  superior  frequency  isolation 
and  stability.  Frequency  isolation  is  rec  .ed  for  laser  radar  designs. 
The  purpose  of  the  work  was  to  produce  a  user  friendly  computer  model 
which  provides  a  realistic  simulation  of  an  injection  seeded  pulsed  C02 
laser.  The  computer  language  chosen  was  QuickBASIC®  to  allow  easy 
future  modification.  QuickBASIC  cannot  compute  complex  numbers  in  the 
manner  of  FORTRAN.  This  limitation  is  overcome  by  the  separation  of 
real  and  imaginary  parts  of  wave  equations.  The  program  will  develop 
results  for  a  set  of  physical  conditions.  The  model  will  not  develop 
the  set  of  input  conditions  required  to  produce  a  specific  result. 

Scope  of  the  Work 

The  computer  program  requires  less  than  128  kilobytes  to  run, 
utilizes  an  EGA  display,  and  does  not  require  extended  or  expanded 
memory.  The  model  is  based  on  a  point  model  of  the  laser.  This  model 
utilizes  four  rate  equations  for  a  four  state  laser.  The  rate  equations 
describe  the  populations  in  the  three  upper  states  and  the  photon 
density  of  the  natural  longitudinal  modes.  The  rate  equations  are 
normalized  to  eliminate  the  dependence  of  the  rate  equations  on  physical 
units.  In  addition,  two  equations  describe  the  cavity  detuning  angle 
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and  power  evolution  for  an  injected  laser.  The  model  ignores  higher 
quantum  levels  and  focuses  on  the  P(20)  transition.  The  number  of 
longitudinal  modes  is  limited  to  a  frequency  range  of  eight  free 
spectral  ranges  higher  or  lower  than  the  transition  frequency.  This  is 
due  to  the  size  of  video  display.  The  video  display  depicts  the 
evolution  of  the  power,  gain,  and  energy  of  any  mode  in  comparison  to 
either  the  most  powerful  mode  or  to  the  summed  power  and  energy  in  the 
laser  pulse.  The  injected  mode  display  shows  the  injected  phase  as 
well. 

Chapter  II  describes  the  model  and  basic  rate  equations.  Chapter 
III  describes  the  modifications  of  the  basic  model  for  multimode 
operations.  Chapter  IV  details  the  considerations  used  in  developing 
the  computer  model.  Chapter  V  contains  the  literature  review  relevant 
to  laser  injection  locking  and  Chapter  VI  is  a  discussion  of  the 
results . 
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inapter  11  aquations  ana  rcoaei 
This  chapter  explores  the  physical  model  on  which  the  computer 
code  is  based.  As  part  of  this  model  the  discussion  will  include  a 
short  description  of  the  spontaneous  emission  term,  Ws.  The  transition 
rates  between  states  and  the  pump  rates  for  those  states  are  delineated. 
Cavity  lifetime  is  defined.  The  injection  field  equation  is  presented. 
The  expression  for  instantaneous  phase  and  injected  signal  field  change 
is  derived  in  a  manner  similar  to  other  works.  A  comparison  between 
steady  state  flux  and  field  equations  verifies  the  validity  of  the  mixed 
approach  instead  of  using  field  equations  exclusively  as  Lachambre  did 
in  his  works  (5:757). 

Model  of  the  Physical  Process 


A 


The  C02  laser  process  involved  is  the  exchange  of  energy  between 
the  asymmetric  stretch  level  of  00°1  and  the  symmetric  stretch  mode 
10°0.  Gilbert  et  al  (3:2524)  proposed  a  model  of  this  process  using  the 
model  shown  as  Figure  1.  Figure  1  shows  the  existence  of  both  the 
higher  energy  and  lower  energy  states  of  C02  as  well  as  the  role  of 
nitrogen  as  a  reservoir  of  energy  for  the  lasing  process.  The 
transition  between  nitrogen  and  carbon  dioxide  is  endothermic.  The 
higher  carbon  dioxide  state  is  actually  more  energetic  than  that  of  the 
nitrogen  by  18  cm-1  .  This  amount  is  small  in  comparison  to  room 
temperature  energy  -  208.3  cm'1  (12:261).  The  population  rate  equations 
that  are  applicable  to  this  model  are  shown  in  equation  set  (1)  (9:4). 

dN. 

=  oc,MNb-Na>  -Y.N.  +  Ycco,Nc-Yci.,N.+Wa 

=  aci|r(Na-Nb)  +  ^Na-YbNb+Wb  (1) 

dN„ 

-  Y cH; ^a- Y cC02 Nc+Wc 

The  symbols  follow  the  international  system  for  units  (S.I.  units) 
unless  otherwise  noted:  c  is  the  speed  of  light,  a  is  the  stimulated 
emission  cross  sectional  area  of  the  C02  molecule,  i|»  is  the  photon  flux 
density  in  units,  Wa,  Wb,  and  Wc  are  the  pumping  rates  for  each 
population,  and  Na,  Nb,  and  Nc  represent  the  concentrations  of  those 
populations.  The  y  symbols  represent  the  collisional  relaxation  rates 
for  the  transitions  between  the  states  noted  in  Figure  1 .  The  reduction 
by  1/4  in  the  equation  for  Nb  is  due  to  the  fact  that  the  transition 
from  the  higher  to  the  lower  energy  state  has  a  four  different  possible 
changes  in  state.  Since  we  are  interested  in  only  one  of  these  changes 
we  incorporate  this  diversity  by  reducing  the  second  term  by  1/4.  In 
addition,  the  rate  equation  for  the  photon  density  is  shown  in  Eq  (2): 


=  4»ac  (Nd-Nb)  _-p^— 'l,+NaWi 


^CAV 


(2) 


where  tCAV  is  the  cavity  lifetime  and  W3  is  the  term  we  will  use  to 
describe  the  spontaneous  generation  of  flux  in  the  manner  of  Gilbert 
(3:2525). 


Cayity.  .Lifetime... -—Lcav 

For  the  purposes  of  this  work  the  cavity  lifetime  will  be  as  shown  in  Eq 
(3)  where  Lres  is  the  length  of  the  cavity  and  R  is  the  reflectivity  of 
the  output  mirror  with  the  assumption  that  the  other  mirror  is  lossless. 


tcAV  = 


2Lres  1 

0  Hi) 


(3) 


In  some  works  (notably  Lachambre),  the  cavity  lifetime  includes  a  factor 
for  transmission  losses  or  absorption  losses  inside  the  cavity  itself. 
Lachambre  defined  the  cavity  lifetime  as  tCAV  =  2r/[ln  ( 1 / ( RT ) ]  where  T 
is  the  transmissivity  of  the  empty  cavity  and  2r  is  the  round  trip  time 
of  a  photon.  T  accounts  for  the  internal  losses  of  the  cavity.  Cavity 
losses  affect  the  photon  density  on  every  trip  through  the  cavity.  An 
equivalent  factor  is  possible.  The  photon  density  is  reduced  by  an 
equivalent  percentage  for  each  time  interval.  This  model  will 
incorporate  this  concept  in  Chapter  III. 

a  -  Effective  Radiative  Cross  Sectional  Area 

The  effect  of  pressure  broadening  is  to  reduce  the  effective  radiative 
cross  sectional  area.  This  reduction  also  depends  on  the  spectral 
separation  from  the  laser  transition.  The  effective  radiative  cross 
section  of  the  molecule  is  Eq  (4)  (12:147;  13:57): 
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(4) 


X2 

®=gfs  (v0,v)  Aco,££ 

where  40  is  the  transition  wavelength.  The  Einstein  A  coefficient, 
ACoeff>  depends  on  the  isotopes  involved  and  is  set  at  either  0.174  for 
C12  or  0.208  for  C13  (4:534;  7:1091).  The  frequency  for  the  transition 
has  an  isotope  dependence;  28.306  THz  for  C12  and  26.889  Thz  for  C13 
(13:23-49).  The  factor  S(v0,  v)  designates  the  line  shape  and  depends 
on  frequency  of  the  transition,  v0,  and  the  frequency  of  interest,  v. 
This  factor  is  normalized  according  to  Witteman  (13:59)  and  Verdeyen 
(12:142)  as  Eq  (5): 


f  S(v0,v)dv  =  l  (5) 

0 

The  line  shape  is  Lorentzian  for  high  pressures.  This  results  in  Eq  (6) 
(13:60)  for  line  shape: 


S  (v0< 


(6) 


The  pressure  broadening  factor  is  taken  from  Witteman  (13:61)  shown  in 
the  next  equation  as: 


aVp-7.58  (  5co,  + ^  5Nj+ .  64  5He+ .  38  $  Hj0+ .  32  $H2  ) 


x  P  •  760 


<  in  MHz  ) 


(7) 


where  £C02,  £N2,  £He,  £H20,  and  £K2  refer  to  the  fraction  of  the  gas  mix  of 
that  component  (in  total  of  1),  the  pressure,  P,  is  in  units  of 
atmospheres,  and  the  temperature,  T,  is  in  Kelvin. 


Transition  Rates  between  States 

The  rates  for  the  various  transitions  are  dependent  on  the  isotope 
of  the  carbon  in  the  molecule,  the  gas  pressure  and  the  relative  amounts 
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of  each  gas  component.  For  C12  the  set  of  rates  from  Lachambre  (5:756) 
using  the  same  notation  as  the  previous  equation  is  as  follows  (9:4): 

Ya=  P  •  760  •  (350  SCo2  +  106  +  85  5„9  +  24000^  +  3864  5H>  > 

YcN,=  P  •  760  •  19100^,  () 

Yocor  P  •  760  •  173705CO2 

Yb=  P  •  760  •  (194  5COj  +  6 5 0  +  32705He  +  4500005HjO  +  710005Hj  ) 

The  rates  for  the  transitions  for  C13  are  shown  in  Eq  (9): 

Ya=  P  -760  •  (1000  5COj  +  304  +  90$He  +  24000^  +  3864  ) 

YC(I,=  P  •  760  •  17370§Ni  (  ) 

Yccor  P  -760  •  191005COj 

Yb=  P  -760  •  (194 5co.  +  650 +  32705He  +  4500005HjO  +  710005Hj  ) 

The  rates  for  the  transitions  for  the  lower  state,  Yb>  have  been 
measured  by  an  induced  fluorescence  technique  where  the  molecules  of  a 
passive  cell  are  excited  and  then  spontaneously  decay  (13:76). 

W-  -  Spontaneous  Generation  of  Flux  Density 

One  can  incorporate  the  effects  of  the  Boltzmann  distribution. 

The  effective  stimulated  emission  cross  sectional  area  is  a  product  of 
the  stimulated  emission  cross  sectional  area  and  fu,  the  thermal 
equilibrium  probability  of  occupation  (3:2535).  The  value  of  fu  is 
roughly  7%  for  the  P(18)  to  P(22)  transitions  at  standard  temperature 
and  pressure.  For  the  P(19)  transition  the  value  of  this  reduction 
expressed  as  fu  is  0.0715  (3:2535).  The  spontaneous  generation  of  flux 
depends  on  the  product  of  the  fraction  of  the  molecules  in  the  upper 
state,  Na,  and  Ws.  Ws  is  a  product  of  the  inverse  radiative  lifetime  of 
the  lasing  transition  (the  Einstein  A  coefficient)  and  the  fraction  of 
photons  that  are  radiated  into  the  small  aperture  of  the  secondary 
mirror  within  the  spectral  width  of  the  axial  mode,  F  (3:2523).  The 
quantity  F  is  shown  in  Eq  (10)  in  derivation: 
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=  s  (v)  civ  Acoe££ 

Xl  fu  avp _ i 


Area  2n 


£u  AVp 


(v-v0 )2+i 


av. 


dv  A, 


coeff 


(10) 


Area  2n 


(v-v0 )2+ 


'iV 

2 


2  2  n 


=  4- 


Area  tCAV 


where  Area  is  the  area  of  the  secondary  mirror,  dv  =1/(2 jr  tCAV) ,  and 
Avp  is  substituted  from  Eq  (7). 


Pump  Rates  for  Each. .State  and  Pump  Cycle 

The  model  of  this  laser  is  not  dependent  on  the  pumping  mechanism. 
This  model  does  not  incorporate  any  accounting  for  penetration  depth  or 
nonhomogeneous  distribution  of  energies  throughout  tne  medium.  Using 
the  ideal  gas  law  the  pumping  rate  is  calculated  *0  or.d  in  2q  (11) 
as : 


W0  =  Effpump  •  9 .65  •  102‘  •  P 

Wa  =  j  Effpump  5COj  -9.65-  102<  •  P 

J  1  vpuinp 

Wb  =  wa 


(11) 


where  the  efficiency  of  the  pumping  mechanism  is  expressed  as  Effpump, 
and  Tpump  is  the  length  of  the  pump  pulse.  These  equations  incorporate 
the  assumption  that  excitation  of  nitrogen  is  3  times  as  effective  as 
C02  and  that  the  pumping  is  indiscriminate  among  the  states  (both  Na  and 
Nb  are  equally  pumped). 


Degeneracy  Ratio 

The  degeneracy  ratio  effects  the  small  signal  gain.  This  study 
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a  rotational  change  of  -1  (13:19).  This  rotational  change  from  the 


upper  state,  j',  to  the  lower  state,  j,  generates  a  degeneracy  as 
(2j'+l)/(2j+l)  =  gu/gi  =  39/41  =  0.95122  (13:63).  Degeneracy  impacts 
the  stimulated  transition  term  in  Eq  (2)  as  a  factor  reducing  the  amount 
of  population  in  the  upper  state  available  to  support  the  transition. 

Field  Equations  for  Injected  Signals 

For  an  amplified  electromagnetic  wave  propagating  in  a  cavity  to 
have  consistent  boundary  conditions  the  length  of  the  resonator  should 
be  n(2/2)  where  n  is  any  positive  even  integer.  The  wave  builds  without 
degradation  due  to  inconsistent  boundary  conditions.  Each  wave  exists 
on  what  we  shall  call  a  natural  longitudinal  mode.  For  an  injected 
signal  not  on  these  resonance  frequencies  a  phasor  analysis  is  done. 

For  such  a  field  the  signal  at  the  end  of  a  round  trip  (designated  2  r 
in  time)  we  have: 


-  -  (12) 
Ei(t  +  2t)=  gainroundtrip Et( t )  +E0(t  +  2t) 

where  E*  is  the  evolving  circulating  signal  and  E0  is  the  injected 
signal.  The  round  trip  power  gain  of  a  cavity  is  G  given  by  Eq  (13): 

G  =  Tes(c)1  =  g2  (13) 

where  T  is  the  cumulative  transmission  gain  (0<T<1)  through  the  passive 
cavity  for  a  round  trip  not  including  output  coupling,  a(t)  is  the  power 
gain  per  unit 'length  of  amplifier  and  1  is  the  length  of  the  active 
section  inside  a  cavity  of  length  Lres  becomes 

gain,roundcriF)  =  gRe2jVf—  (14) 

After  the  incorporation  of  Eqs  (13)  and  (14),  the  round  trip  gain  for 
any  injected  signal  may  be  expressed  as  (10:395): 
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9ain  (round  crip)  (RT)2expl 


(15) 


where  R  is  the  product  of  the  reflectivities  of  the  mirrors  in  the 
optical  path.  The  offset  or  detuning  angle  is  the  difference  between  a 
nearest  natural  mode  and  the  injected  signal.  This  is  the  difference 
between  2kLres  where  k  is  the  wave  number  of  the  injected  signal  and  2mr 
where  n  is  an  integer  large  enough  to  make  the  difference  to  a  range  of 
values  between  -it  and  it.  We  will  designate  this  difference  as  dd 
(10:387).  We  convert  the  effects  of  the  reflectivity  and  transmission 
losses  from  Eq  (13).  Next  we  take  the  natural  logarithm  to  each  side  of 
Eq  (15).  We  note  that  2r  is  equal  to  2  *  Lres/c.  After  dividing  by  2r 
an  expression  for  the  exponential  representation  of  the  evolving  field 
is  Eq  ( 16) : 


In  (gainroundtrip)  _  (  -1  c  ct  ( t )  l-cd8^ 

2 1  '  (2t„v  4  L  J  2L  )  (16) 

=  m  ( t ) 

A  differential  equation  for  the  behavior  of  this  wave  is  then  developed 
through  a  Taylor's  expansion.  The  Taylor's  expansion  for  ex  where  0<x<l 
as  ex  =  1  +  x  +  x2/ 2!  +  x3/3!  +  ...  When  x  is  small  the  first  two  terms 
of  the  expansion  are  a  sufficient  approximation.  For  a  small  gain 
scenario  the  value  of  a  circulating  wave  at  the  end  of  a  round  trip  is 
(10:395-396): 


-  2t  -S;  _ 

Et(t+2t)=e  dtEt(t) 

By  defining  grt  as  the  round  trip  gain  in  Eq  (18) 

a  =«2Tm(c) 

arc  w 


(17) 


(18) 


we  replace  the  round  trip  gain  in  Eq  (12)  with  Eq  (18): 


(19) 


e^Ei  (  t)=e2tra(t)Ei  (t)+E0  (t  +  2t) 

Through  manipulation  and  using  the  forementioned  Taylor's  expansion  as 
an  approximation  we  develop  a  differential  equation  for  an  evolving 
field  as: 


a 

dt 


E*  (t) 


m  ( t )  E*  ( t ) 


f  m  ( t )  E0  ( t »  2 1 ) 


(20) 


Comparison  of  Flux  versus  Field  Equations  _ in  Steady..  State 

The  flux  density  rate  equation  is  equivalent  to  the  field  equation 
in  the  steady  state.  The  flux  rate  Eq  (2)  is  now  written  as: 

4^  =  qr^-t+osc  (AN) i|»  (21) 

CCAV 


where  the  degeneracy  factor  gu/g1  is  incorporated  as  seen  in  Eq  (22): 

AN  =  Nd  -  —  •  Nb  (22) 

From  the  field  model  just  derived 


4§-m(t)E(t)  +m(t)  -Jggg— 

dt  e  ( c  J  _i 


(23) 


where  E  is  the  evolving  field  of  the  longitudinal  mode  and  Espont  is  a 
terra  for  the  spontaneous  generation  of  the  field.  After  subtracting  the 
spontaneous  terms  the  field  equation  is 

■jjf  =m(t)E(t)  (24) 

dt 

The  steady  state  for  the  equation  of  m(t)  for  natural  longitudinal  modes 
with  no  offset  angle  m(t)  is  Eq  (25): 
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(25) 


m(C)=_zi_  +  .£l£ili. 

2  4  L 


With  the  power  gain  coefficient,  a(t),  as  follows 


«(t)  =  a9  (H.  -  I*  *Nb)  •  Q 


' 9  '“d  77"  "b  '  i 
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(26) 


m(t)  becomes 


m  ( t) 


o8ANc  i 
2  2  t.-^v 


(27) 


The  exponent  in  the  field  equation  is  dimensionless  and  any  quantity  may 
be  used  to  define  the  field  E.  Therefore,  defining  the  field  as  the 
square  root  of  the  flux  density  the  derivative  of  the  field  is: 

_i 

dE_  dj;  (28) 

dt~  dt  ~2V  dt 

Eqs  (24)  and  (25)  with  Eq  (28)  simplify  to: 

i  (29) 

CCAV 

Thus,  the  two  forms  for  photon  density  and  field  (21)  and  (24)  are 
equivalent  as  expected.  An  injected  signal  requires  a  phasor  analysis. 

Injected  Field  and  Phase  Equations 

This  section  relies  on  the  derivation  of  two  first  order  non¬ 
linear  differential  equations  as  published  in  Computer  Modeling  of  Gas 
Lasers  (10:395-400): 

E0=E0(c)ej“^  (30) 


Let  the  value  of  the  injected  signal  be  represented  by  the  expression 
above .  The  previous  srguinent  shoved  in  the  slowly  vsrying  ene lysis  when 
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the  injected  signal  was  not  changing  quickly  in  time,  i.e.  E0(t) 
E0(t+2r); 


dEi  m  ( t )  E0  ( t ) 

— i=m(t)  Ei  (t)  + - 5 


dt 


e  2tm  ( t) 


(31) 


The  equation  for  the  evolving  injected  field  is  divided  into  magnitude 
and  phase: 


El(t)=|E1(t)|e^9(t> 


(32) 


(32)  becomes  after  differentiating: 


dt 


dt 


(33) 


This  in  turn  is  the  instantaneous  field  represented  below  (10:395): 


|E,  (t ) 


dt 


(34) 


The  value  of  m(t)  in  Eq  (16)  is  treated  by  dividing  the  expression  into 
real  and  imaginary  parts  so  that  m(t)  -  x  +  jy.  The  expression  for  m(t) 
is  substituted  into  Eq  (31).  Euler's  identity  is  used  on  the 
exponential  term  in  Eq  (31).  With  the  additional  notation  for  the  new 
variables  c  and  d  seen  in  Eq  (35) 


-1  ,  eg  ( t )  1 


2  t 


CA7 


4L 


cd8 

2L 


(35) 


c«=e2txcos  (2xy)  -1 
d=e2tx sin  (2ty) 


we  can  compare  like  terms  and  derive  the  same  equations  in  literature 
for  phase  and  magnitude  (10:400). 


(36) 


dri  _  +  cos 8  (xc <-yd)  -sin8(xd-yc) 

dt  ~  ^  cz+d2 


d8  _  1  sin0(xc+yd)  -cos0(xd-yc) 

dt~y~n  c2+d2  u/; 

The  combination  of  Eqs  (36)  a’ d  (37)  describes  the  evolution  of  an 
injected  signal  and  like  the  population  and  flux  equations  previously 
detailed  are  coupled  differential  equations.  However,  Eqs  (36)  and  (37) 
are  nonlinear  differential  equations  due  to  the  transcendental  terms  for 

e. 


Runge-Kutta  Integration  and  Multimode  Operations 

The  fourth  order  Runge-Kutta  method  of  integration  is  accurate  to 
the  fifth  order  of  the  binomial  expansion.  Numerical  integration 
provides  approximate  solutions  to  differential  equations.  For  all  the 
populations  and  flux  densities  given  t0  as  a  start  of  an  interval  of 
duration  h: 


Na(t0  +  h)  «  Na(t0)  +  J  f  (Na,Nb,Nc,  flux,  time)  dt  (38) 

c0 

The  end  value  of  any  population  or  flux  equation  depends  on  the 
initial  value  of  that  variable  and  the  effects  of  change  induced  on  all 
of  the  variables  during  the  interval.  Eq  (38)  is  applicable  to  the 
other  populations  and  the  flux  density  as  well.  A  series  of  successive 
approximations  in  Taylor  Series  expansion  becomes  equation  set  (39)  for 
a  generic  variable  y  that  is  a  function  of  x  and  y: 
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(39) 


y  i  t0+h)  =  y  (t0)  +k 

where  k  is  found  by  approximation  as  follows: 
ki  =  hfQ 

k2  =  hf  (x0  +  -|h,  Yo  +  ’lki) 

k3  »  hf  (x0+-|h,  y0  +  -|k2) 
k4  =  hf  (x0+h,  y0+k3) 
k  -  i  (k1  +  2k3  +  2k3+k4) 

O  , 

The  values  of  the  integrand  are  relatively  insensitive  to  changes 
in  other  variables  in  small  gain,  small  interval  models.  However,  there 
are  many  potential  longitudinal  modes.  If  each  mode  is  not  selected  in 
isolation  from  other  modes  we  have  designated  an  order  or  precedence. 
This  is  not  the  case  in  reality.  Instead  we  must  hold  the  populations 
constant  for  all  the  flux  calculations.  Each  mode's  flux  builds  in 
isolation  from  other  modes,  affected  by  the  "instantaneous"  population 
only.  These  individual  fluxes  are  then  summed  to  obtain  a  new  total 
flux  for  population  effects.  At  the  next  time  step  the  integration  of 
individual  mode's  flux  is  based  on  the  previous  flux  of  that  mode  and 
again  summed.  This  process  is  shown  in  Figure  2.  Thus,  the  flux  for 
each  mode  is  based  on  the  instantaneous  populations  frozen  in  time. 

This  method  cf  integration  is  accurate  to  the  fifth  order  of  the  power 
series  expansion  of  the  function  used  to  describe  the  independent 
variable  (1:492). 


Figure  2  Computation  Algorithm  for  Runge-Kutta  Approximation 
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This  chapter  explains  the  extent  to  which  the  basic  equations  are 
modified  for  computer  modeling.  The  first  of  these  modifications  is  the 
normalization  of  the  flux  and  population  equations  to  dimensionless 
units.  The  second  modification  is  the  inclusion  of  the  line  shape 
factor  in  the  flux  equations  to  account  for  reduced  cross  sectional  area 
due  to  frequency  separation  of^the  various  longitudinal  modes  from  the 
transition  frequency.  The  assumptions  inherent  to  the  model  and  the 
simulation  default  settings  are  -discussed. 

Normalization 

Normalization  is  a  preferred  technique  that  eliminates  physical 

units  from  an  equation  or  a  system  of  equations.  Additionally, 

normalization  has  the  added  benefit  of  eliminating  computational 

problems  involved  in  the  Runge-Kutta  technique.  Without  normalization, 

•** 

the  chosen  integration  technique  suffers  a  major  error  at  the  start  of 
computations.  The  various  populations  are  low  in  comparison  to  the 
decay  rates  and  a  negative  population  density  develops.  Although  this 
effect  can  be  countered  by  resetting  the  populations  to  zero  density  it 
is  more  practical  to  normalize.  The  electromagnetic  field  equation  may 
be  in  any  units  desired.  The  equivalent  expression  for  power  in  the 
photon  rate  equation  is  in  units  of  photon  density  per  second.  The 
electromagnetic  field  is  directly  related  to  the  square  root  of  the 
power  produced.  As  the  power  produced  is  directly  related  to  photon 
density  per  second,  the  electromagnetic  field  equation  defined  in  units 
of  the  square  root  of  photon  density  per  second  is  an  equivalent 
expression  as  shown  in  Chapter  II.  The  equations  are  normalized  by 
choosing  the  critical  photon  density  as  the  number  of  choice.  This 
number,  Dt,  taken  from  Gilbert  (3:2528)  is  Dt=l/(tCav*acenter*c)  in 
meters"3.  The  populations  and  flux  variables  are  divided  by  D..  Next 


the  pumping  rates  are  multiplied  by  cavity  lifetime,  tCAV,  and  divided 
by  Dt.  Finally,  the  transition  rates  are  multiplied  by  tCAV.  This 
results  in  the  series  of  equations  shown  as  equation  set  (40)  (3:2528). 

<&  =  ± 

Dt 

_  _  N(a,b,c) 

n<a.b.c>  5^—  (40) 

r<a,b,c)  =  Y(a,b,c)  ‘tCAV 

-  W  .  tCAV 

Wd,b,c,s  “  w{a,b,c,s) 

The  complete  equation  set  modified  from  Eqs  (1)  and  (2)  in  Chapter  II  is 
shown  as  a  matter  of  record  as  set  (41)  below.  .  ■ 


- j7~  =  *(nb-na)  -rana+rcC0  nc-rcNjna+wa 


=  <D(na-nb)  +-Sna-rbnb+w, 


dna 

dt 
dnb 

dt 

dnc  n  n 

=  ^cHJna_^ocoJnc+wc 

d4*  =  ${na-nb)--^-«+naw. 


(41) 


dt 


^CAV 


Line...S.hap..e^.a.q.t.o.r-in.  J;hs.J.lu^,.,,RatQ..,EquatiQii 

Eqs  (5)  and  (6)  in  Chapter  II  detailed  the  line  shape  factor  for  a 
homogeneously  broadened  condition.  Recall  that  S(v0,  v)  was  found  to 
be: 


AV 

S(v0,v)=  — £ 


2lt  (tv  Y 

(v0-v)2+l^2 


(42) 


and  that  the  effective  radiative  cross  sectional  area  was  found  to  be: 


X2 

<*=8^  S  (vo'v)  Aco *f£ 


(43) 


The  formulation  of  the  "line  center"  stimulated  emission  effective 
cross  sectional  area  is  derived  by  considering  when  v  =  v0  exactly. 

When  this  condition  is  met  S(v0,  v)=  2/in  Avp)  and  so  incorporating  rhe 
distribution  from  Chapter  II,  Eq  (10): 
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^center" 


f  X2 
•LuAo 


8tc  tcAvp  eoe£f 


(44) 


This  gives  in  turn  an  effective  radiative  cross  sectional  area, 
^effective,  as  seen  in  Eq  (45): 


®9ff9ctive“ac9nter*  2  ? (45) 

This  factor,  elective*  wiH  be  the ‘factor,  ae,  used  in  the  flux 
rate  equations.  Recalling  that  Dt  incorporates  <7e,  we  must  modify  the 
cross  sectional  area  in  the  flux  equation  to  account  for  the  frequency 
separation  of  the  longitudinal  modes. 


Absorption  Losses 

The  variable  T  in  Eq  (13)  of  Chapter  II  accounts  for  the 
transmission  losses  on  a  round  trip  basis.  In  accordance  with  the 
desire  to  account  for  "productive"  output  power  the  cavity  lifetime 
described  in  Eq  (3)  of  Chapter  II  does  not  account  for  those 
transmission  losses  in  the  cavity  due  to  absorption.  These  losses  in 
the  cavity  are  accounted  for  by  reducing  the  photon  flux  during  each 
interval  by  an  indexed  amount  for  all  the  natural  longitudinal  modes. 

The  injected  signal  will  be  degraded  by  an  equivalent  amount  with  the 
accounting  for  the  difference  between  field  and  power.  For  the  natural 
modes  during  a  time  period,  t,  for  a  cavity  with  a  round  trip  time  of  2r 
and  a  distributed  loss  photon  attenuation,  0,  the  factor  to  multiply  the 
flux  density  by  is  designated  Lossinterval  as  given  by: 


Loss 


“  1  — 
interval  ~  -1-  2  T  ^ 


(46) 


Lossinterval  accounts  for  any  size  time  period  even  if  the  interval  is 
chosen  to  be  very  small  order.  The  loss  for  the  electromagnetic  field 
of  the  injected  signal  at  each  round  trip  is  simply: 
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l°ssinMrval  =  1  -  #  (47) 

The  interval  size  is  forced  to  that  of  a  round  trip  time  when  an 
injected  field  is  studied  and  an  indexed  expression  for  field  losses  is 
not  needed.  This  avoids  a  complicated  equation  set  that  would  include 
accounting  for  positional  field  evolution  of  the  injected  signal  in  the 
cavity.  Since  this  a  point  model  of  a  laser,  we  will  not  force  this 
derivation. 

Assumptions  Used  With  the  Model 

Spatial  hole  burning  or  inhomogeneity  due  to  physical  location  is 
not  considered.  Another  assumption  is  the  lack  of  transition  between 
nitrogen  and  ground  state  and  the  grouping  of  the  other  rates  to  include 
all  transitions  out  of  that  state.  In  addition,  in  this  point  model  the 
round  trip  gain  must  be  smaller  than  one.  This  is  due  to  the  Taylor's 
expansion  limitation  on  the  argument  described  in  Chapter  II  on  page  11. 
As  Gilbert  (3:2525)  explained  these  are  reasonable  assumptions  if  the 
length  is  short  and  the  other  physical  assumptions  previously  described 
are  valid. 

Default  Conditions 

For  further  study  a  set  of  default  conditions  is  required.  A 
short  description  is  given  below.  In  addition  to  common  physical 
constants  for  the  speed  of  light  and  Planck's  constant,  the  degeneracy 
ratio  and  the  percentage  of  molecules  in  the  P(20)  transition  are  not 
changed  throughout  the  simulation.  The  Einstein  A  coefficient  is  set 
for  C12  at  0.174  sec'1.  The  length  of  the  resonator  is  1  meter.  The 
reflectivity  is  set  to  71.636  percent.  The  area  of  the  secondary  mirror 
is  0.0004  square  meters.  The  gain  section  has  a  gas  mixture  of  10%  C02, 
10%  N2,  and  80%  He  confined  at  1  atmosphere  pressure  at  300°K.  The  pump 
pulse  is  rectangular  in  shape,  has  an  efficiency  of  0.2  and  the  duration 
of  the  pump  is  10  cavity  lifetimes.  The  laser  pulse  output  is  tracked 
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to  400  cavity  lifetimes.  The  cavity  is  not  Q  switched  and  nominal 
losses  are  3%.  During  the  pump  pulse  and  until  30  cavity  lifetimes  the 
computational  interval  size  is  0.1  cavity  lifetimes.  From  30  cavity 
lifetimes  until  400  cavity  lifetimes  the  interval  size  is  1  cavity 
lifetime.  The  normal  number  of  longitudinal  modes  tracked  is  three.  A 
schematic  of  the  nominal  model  is  shown  in  Figure  3. 
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Chapter  IV  Computer  Program  Implementation  Considerations 

This  chapter  describes  the  computer  on  which  the  program  was  run. 
The  software  limitations  are  detailed.  The  display  screens  are 
described  and  examples  provided.  The  program  is  exercised  to  determine 
the  limits  of  a  noninjected  model  in  regards  to  pressure  and 
computational  step  size. 

Machine  Details  and  Considerations 

The  program  was  written  in  QuickBASIC  4.5®  for  implementation  on 
an  IBM-PC®  or  compatible.  The  choice  of  language  was  not  arbitrary  -  a 
previous  effort  tracking  single  mode  operations  was  already  available. 
The  program  compiles  to  a  final  size  of  120  Kbytes.  The  computational 
run  time  is  under  three  minutes  using  the  standard  defaults  on  a  25  M'nz 
80386  based  computer  with  a  math  coprocessor.  The  speed  of  execution 
depends  largely  on  the  number  of  side  modes  tracked,  the  various  step 
sizes,  and  the  length  of  time  to  track  during  the  pulse  process.  A  math 
coprocessor  greatly  speeds  the  time  of  execution. 

Software  Considerations 

The  program  was  written  so  that  the  majority  of  the  variables  were 
expressed  in  single  precision  format.  This  was  based  on  the  memory 
requirements  for  each  type  of  precision  -  double  precision 
representation  requires  4  bytes  to  represent  a  number  while  single 
precision  representation  requires  only  2  bytes.  The  program  maintains 
in  active  memory  at  each  time  interval:  the  population  of  each  state; 
the  power,  energy,  and  gain  of  each  longitudinal  mode;  and  the  power, 
gain,  energy,  and  the  instantaneous  phase  angle  of  the  injected  signal. 
The  simulation  requires  a  significant  amount  of  active  memory  for  16 
modes  (the  maximum  allowed)  and  a  maximum  interval  size  of  400  cavity 
lifetimes . 
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The  computational  routines  are  not  the  largest  part  of  the 
program.  The  computations  are  mostly  contained  in  one  subroutine 
(RungeKutta4  for  integration)  with  two  function  calls  (differential 
equations:  F  for  the  populations  and  aModeflux  for  natural  modes'  flux). 
An  additional  subroutine  is  called  for  the  computation  of  the 
instantaneous  phase  and  for  the  computation  of  the  change  in  the 
injected  field  strength.  The  majority  of  the  program  code  generates  a 
complex  set  of  display  screens  and  options  that  display  the  following 
normalized  to  the  largest  value  of  each  category  during  the  pulse  or 
alternatively  to  the  peak  value  of  all  the  modes  summed:  (1)  relative 
cross  sectional  areas;  (2)  the  value  of  the  injected  signal  cross 
sectional  area  relative  to  the  cross  sectional  areas;  (3)  the  evolution 
of  the  power,  gain,  and  energy  of  any  natural  longitudinal  mode 
(additionally  phase  for  the  injected  signal);  (5)  the  instantaneous 
power,  gain,  and  energy  profile  of  all  modes  at  any  time;  (6)  a 
numerical  display  of  a  mode's  characteristics  with  population  densities 
for  any  period  of  time;  or  (7)  a  numerical  peak  power  and  peak  energy 
analysis.  Gain  curves  on  the  graphs  are  always  based  on  the  highest 
gain  by  any  single  mode.  An  example  of  the  display  is  shown  in  Figure 
4.  Figure  4  is  a  display  of  all  of  the  modes  at  an  instant  of  time.  At 

the  top  left  is  the  time  in  the  pulse  chosen  for  the  display.  The  top 

right  indicates  that  this  particular  graph  is  normalized  to  the  peak 
values  of  the  total  power  and  total  energy.  Each  individual  mode  has  a 
group  of  three  bars  that  chart  the  value  of  the  particular  field  in 
comparison  to  the  maximum  total  peak  power,  the  maximum  total  energy, 
and  the  gain  for  the  mode  in  comparison  to  the  highest  gain  of  any  mode. 

At  the  right  margin  is  a  listing  of  the  time  of  the  maximum  power  of  the 

most  powerful  mode,  the  time  of  maximum  energy  of  the  most  powerful 
mode,  the  time  of  the  peak  gain  of  the  most  powerful  mode,  the  time  of 
the  peak  total  power,  and  the  time  of  the  maximum  energy.  At  the  bottom 
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of  the  display  is  a  listing  of  the  arbitrary  mode  number.  Mode  numbers 
are  assigned  in  order  of  increasing  frequency  separation  from  the 
transition  frequency.  Mode  one  is  always  the  closest  natural 
longitudinal  mode  to  the  transition  frequency.  The  mode  numbers  are 
assigned  alternatively  with  even  numbers  for  progressively  higher 
frequency  natural  longitudinal  modes.  Odd  mode  numbers  are 
progressively  lower  in  frequency.  The  display  reinforces  this  with 
commentary  notations  at  the  bottom  left  and  bottom  right.  Finally,  the 
far  right  bar  set  is  the  total  power  and  total  energy  at  the  same  time 
as  noted  at  the  top  left.  Figure  4  shows  this  multimode  display  for  an 
non-in jected  case.  On  the  actual  display  the  bars  corresponding  to 
power,  gain,  and  energy  are  all  in  color.  All  individual  power  bars  are 
blue,  all  gain  bars  are  red,  and  energy  bars  are  yellow.  The  total  peak 
power  and  the  maximum  energy  are  depicted  in  purple  and  green 
respectively. 


Figure  4  Multimode  Display  for  Nominal  Conditions  at  30  Cavity 
Lifetimes 

Figure  5  shows  the  same  display  at  400  cavity  lifetimes.  Figure  6 


I 


highlights  the  peak  power  distribution  normalized  to  the  total  power  at 
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the  instant  of  peak  power  for  the  nominal  case.  Figure  7  shows  the  same 
distribution  normalized  to  the  largest  peak  power  of  all  of  the  modes. 


Figure  5  Multimode  Display  for  Standard  Conditions  at  400  tCav 
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Figure  6  Display  at  18  tCav  Scaled  against  the  Sum  Total  Peak  Power 
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15  13  11  9  7  5  3  1  2  4  6  8  10  12  14  16 
LOU  FREQ  Return  to  cont  or  A  for  another  tine  HIGH  FREQ 


Figure  7  Same  Conditions  as  Figure  5  Scaled  Against  Mode  1 

The  gain  components  on  the  multimode  and  time  evolution  displays  are 
always  scaled  to  the  largest  gain  of  any  single  mode  as  shown  in  Figure 
8.  The  time  evolution  display  will  show  a  single  mode's  characteristics 
over  the  life  of  the  pulse  as  seen  for  natural  modes  in  Figure  9. 


Figure  8  Multimode  Display  at  Occurrence  of  Peak  Gain 
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Figure  9  Time  Evolution  Display 


Figure  10  Nominal  Analysis  for  Number  of  Modes 

Two  important  decisions  are  the  size  of  the  interval  and  th» 


number  of  modes  to  track.  ae  is  inversely  proportional  to  the  amount  of 
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frequency  separation  from  the  line  center.  As  <?e  becomes  smaller  the 
critical  inversion  density  increases.  Thus,  a  mode  with  extreme 
spectral  separation  might  have  little  effect  on  total  power  or  energy. 
This  details  a  dilemma:  although  more  modes  would  better  reflect  the 
actual  behavior  of  the  laser,  tracking  additional  modes  requires  more 
computing  time  and  memory.  Additionally,  the  length  of  time  that  the 
pulse  evolution  is  tracked  decreases.  Coincidentally,  the  subinterval 
size  used  in  the  Runge-Kutta  integration  increases.  Still  another 
consideration  is  the  fact  that  the  entire  populations  Na  and  Nb  are 
available  for  all  the  modes.  The  center  modes  will  depopulate  the 
inversion  quicker  than  any  outlying  modes.  With  the  entire  population 
inversion  funneling  through  one  mode  the  total  power  is  larger  than  the 
summation  of  power  for  multiple  modes.  The  number  of  modes  that  are 
required  for  the  computations  to  be  within  10%  accuracy  is  the  primary 
consideration.  Figure  10  shows  that  any  number  of  modes  will  match  this 
requirement.  The  graph  shows  the  percentage  difference  between  the 
number  of  modes  for  total  energy  and  peak  total  power  against  the 
condition  of  tracking  one  mode  only.  A  better  criteria  is  the 
comparison  of  the  minimum  number  of  modes  required  to  closely  match  the 
maximum  mode  case.  Figure  11  is  a  comparison  of  the  summed  peak  powers 
and  total  energy  for  the  1,  3,  and  16  mode  case.  The  interval  step  size 
was  chosen  for  the  stable  solutions  under  0.2  cavity  lifetimes  plotted 
against  the  single  mode  case  using  a  step  size  of  0.005  cavity 
lifetimes.  Figure  11  shows  a  clear  degradation  for  increasing  time 
steps.  The  limitation  for  a  consistent  solution  is  0.19  cavity  lifetime 
for  hi.  For  step  sizes  larger  than  0.19  the  results  fluctuate  from 
interval  to  interval.  The  solution  for  three  modes  with  hl=0.19  is 
within  9.8%  of  the  similar  value  for  nl=0.005  and  16  modes.  The 
limitation  of  this  approach  is  the  focus  on  the  peak  power.  Any 
approach  considering  total  energy  alone  results  in  a  step  size  of  almost 
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arbitrary  size  because  the  peak  power  pulse  has  little  effect  on  the 
total  energy,  less  still  on  the  energy  of  any  one  mode. 
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Figure  11  Differences  in  Peak  Power  and  Total  Energy 
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Figure  12  Increase  of  Power  due  to  Pressure 

Increasing  the  pressure  of  the  system  has  two  effects.  First  it 
increases  the  pressure  broadened  bandwidth.  Secondly,  pressure  affect 
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the  gain  of  the  system  and  hence  the  power  output  in  an  almost  linear 
fashion  from  0.1  to  2  atmospheres  as  seen  in  Figure  12.  Increasing  the 
pressure  also  affects  the  allowable  step  size  as  seen  in  Figure  13.  The 
relationship  between  pressure  and  interval  size  is  not  linear.  The 
interval  step  size  must  decrease  to  accurately  record  the  peak  power  as 
pressure  increases.  In  Chapter  VI  the  validation  of  this  model  is  at  1 
atmosphere.  Accordingly  the  step  size  of  hl=0.19  is  within  the  desired 
10%  accuracy. 


Increasing  the  number  of  modes  weakens  the  assumption  that  the  P(20) 
transition  is  isolated  from  the  P(18)  and  P ( 22 )  transitions.  The 
spectral  separation  between  modes  is  150  MHz  with  a  cavity  size  of  1 
meter.  However,  the  separation  between  the  P(20)  and  P(18)  lines  is 
54283.22  MHz  and  the  P ( 20 )  and  P(22)  lines  are  separated  by  53548.92  MHz 
(13:24).  This  large  separation  precludes  influence  by  populations  in 
another  transition  for  the  maximum  number  of  modes  on  each  side  of  the 
transition  frequency.  Increasing  the  gain  does  effect  the  maximum 
allowable  step  size.  Phenomenologically  speaking  the  higher  pressure 
causes  higher  gain  which  in  turn  increases  the  change  in  populations  and 
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flux  density  during  a  set  step  size.  As  seen  in  Figures  10  and  11  the 
difference  between  1  and  16  modes  for  peak  power  is  well  within  10%  for 
step  sizes  up  to  hl=0.2  for  any  of  the  detailed  pressures.  The  maximum 
allowable  step  size  does  vary  with  pressure  as  seen  in  Figure  13. 
Although  hl=0.19  is  a  sufficient  step  size  for  pressure  under  1.2 
atmospheres  a  smaller  step  size  is  required  at  higher  pressures.  A 
single  mode  case  study  always  has  a  larger  power  than  a  multimode  case 
study.  This  is  not  due  to  the  fact  that  in  a  multimode  study  some 
frequencies  represent  a  larger  potential  energy  in  photon  energy  while 
others  are  smaller.  At  these  frequencies  and  cavity  sizes  the 
difference  is  not  that  large  -  the  difference  between  energy  of  a  photon 
between  adjacent  modes  is  less  than  0.001%.  The  mode  associated  with 
the  largest  ae  will  have  the  largest  gain  in  the  noninjected  case.  The 
discrepancy  in  gain  leads  to  a  larger  flux  density  for  the  associated 
mode.  The  larger  flux  density  stimulates  a  larger  percentage  use  of  the 
available  population  of  the  upper  state.  When  the  total  upper  state 
supports  the  mode  with  largest  ae,  the  peak  power  increases.  The  change 
in  flux  density  from  Eq  (2)  is  directly  proportional  to  the  size  of  the 
ae  and  the  flux  density.  The  spontaneous  flux  density  is  also 
proportional  to  the  ae.  The  initial  flux  also  has  a  direct  effect  of 
the  charge  in  flux  density.  Clearly,  the  largest  ae  has  an  advantage  in 
the  nominal  case.  This  is  not  true  in  the  injected  case.  A  large 
injected  signal  flux  may  influence  the  pulse  evolution  more  than  the  ae. 
In  a  continuous,  homogeneously-broadened  laser,  the  gain  is  eventually 
reduced  to  threshold  levels.  The  mode  closest  to  line  center  has  the 
highest  gain.  This  mode  also  has  the  largest  ae. 


32 


Chapter  V  Literature  Review 

In  this  chapter  the  significant  literature  is  reviewed  in 
chronological  order.  Lachambre  used  a  modified  flux-population  rate 
equation  with  a  conversion  from  field  to  flux.  Siegman  and  Tratt  et  al 
both  agree  with  Lachambre' s  approach  and  validate  his  findings  of  the 
injected  mode  selection  range.  Recently  Cassard  and  Lourtioz  used 
Maxwell's  equations  to  predict  conditions  when  two  or  more  signals  would 
exist  simultaneously. 

Lachambre 's  Approach 

Lachambre  developed  the  first  locking  study  of  TEA  C02  using  a 
combined  population  and  field  approach.  He  described  the  injected  field 
in  the  manner  of  Chapter  III.  His  results  showed  three  major  regions:  a 
region  where  the  injected  signal  had  no  effect;  a  region  where  mode 
selection  existed  for  a  wide  variety  of  phases  and  powers;  and  a  region 
where  injection  locking  would  occur  with  very  small  detuning  angles  and 
very  large  injected  powers.  Even  for  low  powers  an  injection  range  of 
almost  0.4 n  in  phase  was  observed.  He  based  his  efforts  on  a  constant 
ae  for  the  C02  molecule  and  for  a  total  of  12  longitudinal  modes.  The 
results  are  instructive  as  depicted  in  Figures  14,  15,  and  16.  Figure 
14  shows  the  relationship  between  the  initial  detuning  angle  or  offset 
and  the  portion  of  the  total  pulse  energy  due  to  the  injected  signal. 

The  vertical  scale  is  the  portion  due  to  the  injected  signal  E£  over  the 
sum  total  of  the  injected  signal  plus  the  energy  in  the  natural  modes 
Es.  The  injected  signal  has  a  strength  of  5  X  10(7_2m)  W/cm2  where  m  is  a 
number  from  2  to  9  as  detailed  on  the  figure.  The  portion  of  the  total 
energy  due  to  the  injected  signal  is  very  high  for  detuning  angles  under 
36°  for  the  injection  powers  shown.  For  any  signal  injected  with  an 
offset  of  0.5n  or  greater  there  is  little  contribution  from  the  evolving 
field  to  the  total  energy. 
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Figure  14  Mode  Selection  Region  (5:761) 

In  Figure  15  the  injected  signal  exhibits  definite  effects  on  the 
evolving  signal.  According  to  Lachambre  the  total  peak  power  of  all  the 
natural  modes  when  the  detuning  angle  is  0.5n  is  45  MW/liter  with  the 
peak  occurring  at  700  nanoseconds  for  a  total  energy  of  4.4  Joules/liter 
(5:759).  A  strong  injected  signal  causes  an  early  reduction  in  the 
population  inversion.  This  is  seen  in  Figure  16  by  the  earlier  peak 
power  at  600  nanoseconds  and  the  reduced  peak  power  of  29  MW/liter.  The 
total  energy  increases  to  4.5  Joules/liter  due  to  the  longer  power 
pulse.  The  lower  charts  in  Figures  15  and  16  depict  the  value  of  2xd/<{> 
on  the  vertical  axis  where  2r  is  the  round  trip  time,  8  is  the 
instantaneous  change  in  the  phase  and  <f>  is  the  initial  detuning  angle  of 
the  injected  signal  form  the  nearest  natural  mode.  If  the  quantity  2x8 
is  equal  to  4>  then  the  evolving  signal  follows  the  detuning  angle.  This 
is  termed  the  mode  locked  condition  and  the  injected  signal  is  said  to 
operate  in  the  mode  selection  region.  If  2x8  is  equal  to  zero  the 
amplified  injected  signal  is  not  locked  to  the  detuning  angle  and  is 
locked  to  the  amplified  injected  signal.  This  is  the  true  injection 


34 


locked  case.  In  the  low  power  case  we  see  a  beat  frequency  oscillation 
in  Figure  15.  This  is  due  to  the  competition  between  the  mode  and  the 
injected  signal  when  the  selected  mode's  power  decays.  High  power 


Figure  15  Mode  Selection  for  Low  Powered  Injection  (5:760) 
injection  causes  an  earlier  reduction  in  the  inversion  before  locking  to 
lock  to  the  injected  signal.  This  suggests  that  injection  locking  is 


possible  only  when  the  gain  is  low  (below  threshold).  Figure  16  shows 
the  signal  evolution  at  very  small  detuning  angles  with  very  high 
injected  powers.  The  injected  power  peak  is  reduced  even  further  to  9 
MW/liter.  The  peak  power  occurs  even  earlier  at  500  nanoseconds.  The 
secondary  peaks  are  less  extreme  and  the  pulse  tail  is  longer.  The 
signal  is  effectively  injection  locked,  and  the  total  energy  in  the 
injected  signal  increases  to  6  Joules/Liter. 


Tratt  CO?  Laser  Injection  Study 

Tratt  completed  an  extensive  review  of  injection  seeding  in  1985. 
This  effort  detailed  the  active  mode  selection  schemes  which  lead  to  the 
selection  of  one  single  longitudinal  mode  (SLM)  (Tratt:235).  The 
features  of  these  schemes  show  a  suppression  of  the  gain  switched  spike 
of  radiated  power  and  a  longer  tail  over  the  time  of  the  pulse  as 
depicted  in  Figure  17  (11:237). 


As  seen  in  Figure  17  the  introduction  of  a  signal  causes  early 
depletion  of  the  inversion  level  and  leads  to  earlier  signals  and  lower 


powers.  Eq  (48)  is  the  v  -fe  used  in  his  analysis  and  it  similar  to  the 
description  in  Chapter  II. 


E(t)  =  /R  E(t-2T)eC«('>‘^]+iinj80„d 


(48) 


Here  G(t)  is  the  real  part  of  m(t)  expressed  in  Eq  (16)  of  Chapter  II,  R 
is  the  reflectivity  of  the  secondary  mirror  and  4>  is  the  imaginary  part 
of  m(t)  from  Eqs  (12)  and  (16)  earlier.  The  steady  state  solution 
(11:245)  of  Eq  (48)  is  when  G(t)=0.  Eq  (48)  is  separated  into  real  and 
imaginary  parts.  This  results  in  the  real  part  of  the  field,  9),  and  the 
imaginary  part,  3,  as: 


St  [E  (0  )]  -  Einiecced 


1  -\/RCQSC> 
l  +  R-2  v/Rcos<j> 


(49) 


»  [E(0)]  =  Eindected 


v/Rsin4> 


(50) 


l+R-2  v/R  cos<J> 

From  these  two  expressions  Tratt  developed  an  expression  for  the 
round  trip  phase  shift  d(t)  as: 


tan  [0  ( t)  -  8  (t-t )  ] 


g[B(t)]  ft[E(t-t)j  -  *[E(t)]  g[E(t-t)]  (51) 
*[E(t)J  St  [E  (t-x )  ]  +  g  [E  (t)  ]  g[E(t-T>] 


and  the  instantaneous  frequency  F(t)  becomes  (0(t)-#(t-r)  )/<£  (11:245). 
The  implications  are  that  a  zero  phase  shift  indicates  a  locking  to  the 
injected  signal  but  that  a  phase  shift  change  implied  a  resonance 
condition  with  the  nearest  cavity  mode  like  Lachambre.  Tratt  documented 
the  three  regimes  noticed  by  Lachambre  (11:246).  Tratt  noticed  the 
existence  of  the  beat  frequency  oscillations  in  the  higher  power  cases 
due  to  interactions  between  the  existing  cavity  signal  and  the  injected 
signal  (11:246-247). 

Tratt  found  that  injected  signals  with  a  detuning  angle  from  zero 
to  n  had  an  effect  on  the  power  output  profile.  The  effects  became  less 
noticeable  for  larger  and  larger  offsets  as  depicted  in  Figure  18.  Many 
authors  agree  with  the  limitation  of  the  detuning  angle  of  O.^tt  for 
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Figure  18  Injected  Regime  (11:248)  <f>  is  the  detuning  angle 


injection  locking.  However,  other  authors  using  large  cavities  found 
that  injection  locking  could  be  forced  over  the  entire  spectrum 
(11:260). 


Siegman's  Derivation 

Siegman  objects  to  the  practice  of  using  the  term  "injection 
locking"  because  the  injected  signal  provides  a  set  of  initial 
conditions  from  which  the  wave  grows  with  little  regard  to  the  injected 
signal.  He  would  rather  the  term  "injection  seeding"  apply  in  this 
scenario.  In  a  laser  cavity  the  wave  circulating  at  a  time  t,  Ecirc(t), 
would  be  the  sum  of  the  previous  circulating  wave  that  is  just  reflected 


off  the  mirror,  Ej^t),  and  the  injected  wave,  E0(t),  using  the 
expressions  below  (8:1155): 

Eoiro(t)=Ei(t)+E0(t)  (52) 

Ei(t)  =Ecirc  (t  -2t)e  (53) 

where  g  refers  to  the  cumulative  gains  inside  the  cavity  of  the  previous 
wave  ac  time  t-2r,  wosc  is  the  frequency  of  the  nearest  longitudinal 
mode,  and  o)inj  is  the  frequency  of  the  injected  signal.  When  the  cavity 
is  below  threshold  and  the  signals  are  on  resonance  the  vectorially 
added  signal  results  in  a  steady  state  nominal  wave.  The  injected 
signal  makes  up  the  vectorial  difference  between  the  circulating  wave 
and  the  longitudinal  mode.  For  injected  signals  of  low  power  there  are 
distinct  possibilities  for  injection  seeding.  After  separating  the 
variables  and  using  Euler's  identities  for  exponents  (8:1156): 


oirc  [l-e-3cos(2x  (wSnj-u0SC))]+j[e-9sin(2T  («jnj-wesc))]  1 
This  expression  shows  that  for  the  signal  to  build  the  subexpression 
I  2r(winj  -  wosc)l  must  be  less  than  n/2.  When  the  instantaneous  phase 
4>(t)  is  designated  as  Siegman  did  (8:1160): 


(54) 


(j)(t)-  -(<*>inj-0>osc)  T  •— 

then  the  actual  frequency  wi(t)  is  simply: 


(55) 


i*\  _  d<J>(t)  ^ 

Mi  C - dt-  Wlnj 

”  ““inj  +  wosc  +  uinj 


(56) 


From  Eq  (55)  and  using  a  expansion  of  the  equation  in  the  form  of 


f  (x  +  h)  =  f  (x)  +  hf'(x)  +  yj-  f''(x)  (etc)  (57) 

based  on  the  assumption  that  <oi-w0  <<  l/(2r)  Eq  (52)  becomes: 
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2  T 


1  E, 


(58) 


The  implication  from  Eq  (56)  is  that  the  frequency  of  a  mode 
selected  laser  is  almost  entirely  on  the  resonant  mode  wavelength  and 
for  all  practical  purposes  the  injected  contribution  is  buried.  This  is 
phenomenologically  correct  when  the  circulating  wave  is  much  bigger  than 
the  injected  signal  in  a  mode  selected  regime. 


Cassard..  .and  Lourtioz. Approach 

Cassard  and  Lourtioz  approached  the  process  from  Maxwell's 
equations  and  assumed  linearity.  A  wave  equation  of  motion  for  the 
complex  field  E  in  a  medium  with  a  polarization  P  and  current  density  j 
is : 


.rE-2jk(|f  +  ^|E 


2e0c 


E)  =  -j  (l0<i>2  P 


(59) 


where  a  is  the  conductivity  and  co  is  any  arbitrary  mode  frequency 
(2:2322).  After  fixing  the  boundary  conditions  they  developed  a  set  of 
eigenmodes  (2:2322).  Their  solution  provides  for  the  injected  field 
(2:2324)  as  Eq  (60): 

t*2t 

/  P  \  (  C'JWM  /  cr\\ 

E(t  +  2t)  =  |E(t)  +  v'R  e  «  e-2*n# 

where  g (9)  is  the  complex  gain  term  that  includes  resonant  dispersion 
and  6  =  r/n((Oinj-wc)  and  is  equivalent  to  Lachambre's  expression  cd(?/2L. 
In  this  instance  the  injected  signal  is  measured  prior  to  entering  the 
cavity. 

Cassard  and  Lourtioz  objected  (2:2321)  to  the  assumption  that 
Lachambre  made  concerning  the  lack  of  amplification  of  the  injected 
signal.  However,  the  injected  signal  would  only  experience 
amplification  on  one  trip  through  the  cavity.  At  the  end  of  that  trip, 
the  injected  signal  is  added  to  the  circulating  wave.  Since  the 
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circulating  wave  becomes  much  larger  over  time  than  a  low  power 
injection  signal  one  can  discount  the  contribution  of  the  injected 
signal  over  time.  This  validates  Lachambre's  approach.  For  a  high 
power  injection  signal  there  could  exist  an  injected  signal  amplified  by 
one  trip  through  the  cavity  that  adds  vectorially  with  the  evolving 
field.  This  would  cause  the  evolving  field  to  shift  not  only  due  to  the 
offtuning  of  the  cavity  but  also  due  to  vector  addition  with  the 
injected  signal.  There  could  exist  more  than  one  mode  in  the  cavity  but 
this  regime  will  be  confined  to  large  injection  powers. 

Conclusion 

The  combination  flux-population  rate  equation  approach  is  valid 
for  low  power  injection  signals.  The  incorporation  of  line  shape  in  the 
expression  for  cross  section  will  yield  more  accurate  results  than 
Lachambre  for  there  will  be  a  dependence  on  mode  selection  and  spectral 
separation  for  outlying  modes.  Three  modes  is  enough  to  consider  at 
atmospheric  pressures  -  from  the  line  shape  factor  the  relative  size  of 
the  mode  at  8  free  spectral  ranges  for  the  nominal  mix  is  still  70%  of 
the  highest  value  but  has  little  contribution  to  the  total  signal. 


Chapter  IV  Results 

This  chapter  discusses  the  results  of  the  computer  runs.  The 
validation  of  the  basic  model  is  reviewed  and  the  differences  between 
Lachambre's  results  and  the  computer  program  are  explained.  The 
dependence  of  the  evolving  field  on  various  physical  parameters  is 
depicted. 

Validation 

The  simulation  was  run  against  the  original  model  developed  by 
Stone  from  Gilbert's  article  on  C02  Lasers.  There  was  no  difference 
between  the  two  simulations.  Next  the  model  was  run  in  comparison  to 
Lachambre's  results.  The  program  correctly  predicted  the  same  results 
for  the  peak  power  of  45  MW/liter  but  had  a  much  earlier  peak.  There 
are  two  reasons  for  this  discrepancy.  The  first  is  that  Lachambre 
included  in  his  population  equations  a  collision  transition  between  the 
10°0  and  the  02°0  as  a  relaxation  process  for  the  lower  state.  This 
rate  was  dependent  only  on  C02,  N2,  and  He  gas  mixtures.  He  implicitly 
assumed  that  the  relaxation  from  this  state  was  very  fast  and  that  the 
possibility  of  the  reverse  transition  was  not  likely.  Although  a 
relaxation  operation  is  required  from  the  10°0  state  to  keep  a  strong 
inversion  the  method  of  using  a  lumped  relaxation  process  is  more 
accurate.  The  second  difference  is  in  the  generation  of  spontaneous 
flux.  In  Chapter  II  the  derivation  of  the  spontaneous  term  was  based  on 
the  population  of  the  upper  state  and  the  fraction  of  the  spontaneous 
emitted  photons  that  intersect  the  secondary  mirror.  Lachambre  used  a 
different  approach  incorporating  the  concept  of  amplified  spontaneous 
emission  (ASE)  and  determined  a  one-way  photon  noise  density  to  use  as 
his  field  for  the  natural  longitudinal  modes  (5:757).  This  results  in  a 
spontaneous  flux  density  does  not  depend  on  the  population  inversion. 

The  value  for  current  default  parameters  is  109  times  larger  than 
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Lachambre.  Accordingly  our  peak  is  much  sooner.  Since  the  flux  builds 
rapidly  once  established,  there  is  no  effect  on  the  peak  power 
magnitude. 

Additionally,  Lachambre  did  not  include  the  effects  of  a  reduced 
active  element  size  in  his  flux  rate  Eq  (5:758).  He  accounted  for  a 
gain  section  length  inside  a  larger  resonator  cavity  only  in  his 
generation  of  ASE.  He  did  not  account  for  the  possibility  of  the 
smaller  section  length  in  his  field  equations.  The  rate  of  change  of 
the  populations  is  not  dependent  on  the  power,  but  the  density  of  the 
flux  in  the  gain  section  affects  the  change  in  the  upper  and  lower  C02 
states.  This  effect  is  significant.  In  the  manner  of  Milonni  (6:296) 
one  must  add  the  effect  of  the  reduced  active  section  length  in  the 
cavity  to  both  the  flux  equation  for  the  natural  modes  and  the  coupled 
equations  for  the  injected  field.  Figure  19  shows  how  the  peak  output 
power  is  reduced  as  the  length  of  the  cavity  increases  with  a  constant 
gain  section  length. 


Figure  19  is  a  compilation  for  the  default  parameters  with  an 
active  section  length  of  0.5  meters  and  a  reflectivity  of  90%.  The  pump 
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length  and  the  integration  step  size  were  both  adjusted  for  all  events 
so  that  the  pump  duration  was  300  nanoseconds  and  the  integration  step 
size  was  a  constant  3.17  nanoseconds  regardless  of  the  cavity  length. 

In  an  continuous  wave,  non-in jected  scenario,  the  long  term  laser 
frequency  belongs  to  that  of  the  strongest  mode.  In  a  continuous  wave 
laser  under  homogenous  conditions  similar  to  the  present  case,  the  modes 
with  a  smaller  effective  radiative  cross  sectional  area  are  eventually 
quenched  either  by  losses  or  filters  and  only  the  strongest  signal 
survives.  This  is  somewhat  true  in  this  study  as  well.  Five  modes  are 
more  than  adequate  to  describe  the  power  evolution  as  seen  in  Figures  4 
through  7.  In  Figure  5  less  than  3%  of  the  total  energy  comes  from  Mode 
5.  From  a  phenomenological  view  one  might  think  that  smaller  effective 
radiative  cross  sectional  area  would  effect  the  mode  selection  range. 
While  this  is  true  with  small  pressure  broadened  bandwidths  this  is  not 
the  case  at  atmospheric  pressures.  Figure  20  depicts  the  relationship 
between  injection  near  the  strongest  longitudinal  mode  and  the  next 
strongest  mode  for  the  default  parameters. 
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There  is  not  a  significant  relationship  between  the  total  energy 
of  the  injected  signal  for  either  mode  one  or  mode  three  at  atmospheric 
pressure.  This  effect  is  the  same  regardless  whether  injection  is 
higher  or  lower  in  frequency  around  the  chosen  mode  as  long  as  the 
detuning  angle  is  equal.  Figure  20  also  shows  power  and  offset 
dependence  similar  to  Lachambre's  results. 

The  same  effects  that  Lachambre  noticed  for  injected  signals  are 
also  apparent.  There  is  a  distinctive  beat  frequency  oscillation  as  the 
total  power  weakens.  This  effect  is  more  noticeable  for  higher  injected 
powers  and  is  still  more  noticeable  as  the  total  power  weakens.  As  seen 
in  Figure  21,  the  effect  of  the  oscillation  is  low  at  the  beginning. 

This  oscillation  does  not  occur  during  the  peak  power  spike  but  at  some 
later  time.  The  oscillations  become  larger  and  larger.  Lachambre 
noticed  this  as  well  (Figure  15). 
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Figure  21  Display  for  an  Injected  Signal  1  W/cm2 


Figures  22  and  23  show  how  the  phase  oscillations  increase  in 
amplitude  over  time  as  the  injected  signal  becomes  larger.  Figure  22 
shows  that  as  the  injected  signal  is  increased  the  phase  oscillation 
becomes  progressively  larger.  Later  in  the  pulse,  the  effect  of  the 
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Figure  24  Injected  Signal  Near  Mode  4 


.40E+03 


I  Another  noticeable  characteristic  is  that  the  phase  closely 
follows  the  rising  gain  curve  as  seen  in  Figure  25.  Figure  25  shows  the 
evolution  of  an  injected  signal  near  Mode  4  at  a  detuning  angle  of  24 

J  degrees  with  an  injected  power  of  1  Watt/cm2.  The  most  influencing 
actor  in  Eq  (43)  is  the  factor  i).  t)  is  the  ratio  of  the  evolving  field 
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over  the  injected  field.  As  the  evolving  field  becomes  larger,  the 
second  term  in  the  differential  phase  Eq  (37)  becomes  smaller.  This 
results  in  a  dependence  on  the  initial  offtuning  angle  exclusively. 

Figures  26  through  29  show  the  lack  of  dependence  of  the  amplified 


injected  signal  on  the  mode  number  nearest  the  injected  signal. 
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Figure  26  Multimode  Display  for  Injected  Signal  near  Mode  1. 


Figure  26  is  based  on  injection  near  Mode  1  using  the  standard 
conditions  with  a  detuning  angle  of  24  degrees  and  an  injected  signal  of 
1  Watt/cm2.  The  injected  signal  peak  power  clearly  dominates  all 
longitudinal  modes.  Figure  27  is  based  on  injection  near  mode  4  using 
the  same  conditions  as  Figure  26.  The  dominance  of  the  injected  signal 
is  clear  in  Figure  27  as  well.  In  a  mode  selection  regime  all  of  the 
peak  power  is  due  to  the  injected  signal  as  exemplified  in  Figures  26 
and  27. 
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Figure  27  Multimode  Display  for  Injected  Signal  near  Mode  4. 

However,  due  to  the  larger  spontaneous  generation  of  flux  by  natural 
longitudinal  modes  incorporated  in  this  work  the  total  energy 
contribution  by  the  injected  signal  is  not  as  dominating.  This  is  shown 
in  Figures  28  and  29. 
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Figure  28  Display  of  an  Injected  Signal  near  Moce  1  at  4uo  tCAV. 
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For  an  injected  signal  near  Mode  1  the  contribution  by  the  injected 
signal  is  86%  of  the  total  for  the  given  conditions. 
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Figure  29  injection  for  Same  Conditions  as  Figure  27  near  Mode  4 


For  an  injected  signal  near  Mode  4  the  same  injected  signal  strength  is 
only  77%  of  the  total.  The  injected  signal  suffers  from  inconsistent 
boundary  conditions.  A  larger  gain  potential  is  required  for  the 
amplified  injected  signal  than  the  natural  longitudinal  modes.  As  the 
gain  decreases  the  amplified  injected  signal  dies  away.  The  remaining 
population  supports  the  natural  longitudinal  modes.  This  difference  in 
energy  is  reflected  in  Figure  24.  Figures  28  and  29  use  the  same 
conditions  as  Figures  26  and  27. 

Mode  number  has  little  effect  on  the  injection  locking  regime  as 
seen  in  Figures  30  and  31.  In  Figures  30  and  31  we  used  the  default 
parameters  with  an  injected  signal  strength  of  5  KWatts/cm2  at  a 
detuning  angle  of  24  degrees.  Figure  30  is  for  injection  near  Mode  1. 
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Figure  30  Injected  Signal  with  Power  of  5  KW/cm2  near  Mode  1 

Figure  31  is  for  injection  near  Mode  4.  In  both  cases  the  peak  power  is 
at  the  same  time.  The  injection  phase  is  also  the  same. 
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Figure  31  Display  for  Same  Conditions  as  Figure  30  near  Mode  4 


Effect  of  Iniected  Signal  Power  on  Power  Evolution 


An  increased  injected  signal  strength  causes  an  earlier  depletion 
of  the  population  inversion.  Consequently,  as  the  injected  signal 
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strength  increases  the  peak  power  will  decrease  and  the  time  of  the  peak 
is  earlier.  This  is  seen  in  Figures  32  through  34.  Figure  32  shows  the 
power  due  to  the  amplified  injected  signal  for  various  strengths  in  the 
injection  mode  selection  zone  near  mode  one  using  the  standard 


conditions. 


Figure  33  depicts  the  total  power  of  all  modes,  whether 
spontaneously  generated  or  due  to  the  injected  signal  for  a  similar 
interval  as  shown  in  Figure  32.  Notable  is  the  effect  on  the  total 
power  due  to  a  small  injected  signal.  For  an  injected  signal  of  10'6 
W/cm2  the  detuning  angle  is  too  large  to  allow  mode  selection.  However, 
for  a  signal  of  10"*  W/cm2  the  total  power  is  16  MW/liter  (Figure  33) 
while  the  power  due  to  the  injected  signal  is  12  MW/liter.  At  the 
higher  injection  levels  the  power  is  due  almost  entirely  to  the  injected 
signal  with  little  contribution  from  the  other  modes. 
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Figure  34  is  a  close  up  of  mid-pulse  power  evolution  to  show  how  the 
higher  signal  levels  deplete  the  population  inversion  faster  and  cause 
the  power  to  dissipate  faster. 


Figure  34  Total  Power  Evolution  in  Hid  Pulse 
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As  seen  in  Figures  32  through  34  for  very  low  powers  the  effect  on 
the  total  output  is  very  small.  At  higher  powers  the  injected  signal 
has  a  dominant  effect  on  the  power  evolution. 

Behavior  of  Injected  Signal  at  Various  Detuning  Angles 

Lachambre  found  that  the  mode  selection  region  was  as  large  as 
0.471.  This  region  decreased  as  the  injected  signal  strength  decreased. 
He  also  noted  that  the  peak  of  the  injected  signal  was  retarded  in  time 
as  the  offset  increased.  This  model  predicts  similar  effects. 


Figure  35  Peak  Power  as  a  Function  of  Detuning  Angle  in  Degrees 

Figures  35  through  38  are  based  on  the  output  due  to  model  running  with 
Lachambre 's  parameters. 


Table  I  Lachambre' s  Parameters 

Active  Length  =  1  meter  He  =  .7  atmospheres 
L„es  =1.5  meters  C02  =  .15  atmospheres 
T*  =  72%  N2  =  .15  atmosoheres 
R  =  28%  Volume  =  100  cm3 


In  table  I  is  the  list  of  the  applicable  variables  to  the  current  model. 
T  is  the  same  value  in  Eq  (13)  of  Chapter  II.  This  is  equal  to  (1- 
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/3)/100  at  the  end  of  Chapter  III.  R  is  the  reflectivity.  The  volume  is 
translated  into  the  area  of  the  secondary  mirror  multiplied  by  the 
active  section  length.  The  partial  components  for  the  gases  are 
translated  into  their  appropriate  fractions  for  use  in  the  program. 
Figure  35  shows  that  the  total  peak  power  is  very  low  for  small  detuning 
angles  as  much  of  the  energy  is  consumed  in  the  early  degradation  of  the 
population  inversion  by  the  injection  flux  density.  The  power  of  the 
total  combined  flux  densities  reaches  a  plateau  past  50  degrees.  This 
implies  that  the  injected  signal  is  weak.  This  is  indeed  the  case.  The 
effect  of  the  injected  signal  at  large  detuning  angles  is  small.  From 
Figure  13  Lachambre  found  little  total  energy  contribution  past  0.35tt. 


As  seen  in  Figure  36,  the  injected  signal  dies  quickly  outside  the 
mode  selection  range.  The  same  effect  is  seen  in  the  time  of  the  power 
peaks.  The  smaller  the  injection  detuning  angle  is  the  more  effective 
the  injected  signal  becomes  at  inducing  a  response.  The  injected  peak 
power  occurs  earlier  but  the  peak  power  is  lower. 
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Figure  37  Time  of  Total  Peak  Power 


As  the  injected  signal  power  increases  the  peak  power  time 
decreases  corresponding  to  the  early  depopulation  of  the  inversion 
peak  power  occurs  later  as  the  detuning  angle  increases.  This 
corresponds  to  a  decrease  in  the  effect  by  the  injected  signal.  T 
effect  plateaus  in  correspondence  to  the  undriven  case. 
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Figure  3o  Injected  Peak  Power 


Inside  a  limited  range,  the  majority  of  a  pulse's  power  comes  from  the 
injected  signal.  The  proportion  of  the  energy  that  is  due  to  the 
injected  signal  is  significant.  However,  the  amount  of  the  power  due  to 
injection  decreases  quickly  as  the  size  of  the  input  signal  decreases 
and  the  frequency  offset  from  the  nearest  longitudinal  mode  increases. 
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Within  the  confines  of  the  point  model  for  the  laser  and  the  small 
signal  approximation  Lachambre's  equations  predict  the  behavior  of 
injection  mode  locking  in  C02  lasers.  There  are  three  regimes  of 
injection  locking:  a  region  at  large  detuning  angles  where  the  injected 
signal  has  little  to  no  effect  on  the  pulse  forming  process  and 
contributes  little  to  the  total  power;  a  region  where  the  injected 
signal  influences  the  power  evolution  of  the  cavity  and  enhances  the 
output  of  the  nearest  longitudinal  mode;  and  a  region  at  extremely  high 
power  and  small  detuning  angles  where  the  laser  will  oscillate  at  the 
injected  signal.  Indications  of  injection  mode  selection  are  the 
oscillation  of  the  phase  in  synchronous  the  initial  cavity  detuning 
angle.  Symptoms  of  non-mode  selection  are  the  lack  of  contribution  to 
the  power  output  or  early  oscillation  of  the  phase  in  the  injected 
signal  resulting  in  a  beat  frequency  oscillation.  In  the  mode  selection 
regime  phase  tends  to  follow  the  rising  gain.  This  makes  sense  when  Eq 
(37)  is  examined  in  this  new  light.  As  the  variable  n  grows  larger, 
d#/dt  becomes  simply  the  value  of  the  detuning  angle/time.  When  this  is 
integrated  over  time  the  result  is  the  detuning  angle.  Mode  selection 
has  little  dependence  on  the  nearest  mode.  The  primary  factors  in 
determining  mode  selection  are  the  detuning  angle  and  the  injected 
signal  power  in  that  order.  Mode  number  is  a  distant  third.  Increasing 
the  pressure  will  increase  the  gain  and  hence  the  output  of  the  laser. 
Incorporating  a  reduced  cavity  length  decreases  the  output  power  in  an 
inverse  manner  -  as  the  active  section  to  cavity  length  ratio  decreases 
the  power  output  of  the  laser  decreases. 
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A  homogeneously  broadened  C02  laser  is  examined  as  a  point  model 
with  injection  locking.  Three  distinct  behavior  regions  are  discovered 
in  confirmation  of  previous  works  by  Lac'nambre,  Tratt,  and  Cassard.  Low 
power  injection  at  medium-to-low  off tuning  angles  causes  the  nearest 
longitudinal  mode  to  dominate  all  other  modes.  High  power  injection  at 
very  low  offtuning  angles  causes  the  injection  signal  to  amplifiy  and 
dominate  all  of  the  longitudinal  modes.  The  injected  signal  has  little 
effect  on  the  power  or  energy  evolution  of  the  multimode  cavity  at  large 
offtuning  angles.  A  computer  model  written  in  QuickBASIC®  for  use  on  an 
IBM-PC/AT®  or  compatible  is  used  to  examine  the  behavior  of  injected 
signals.  Field  and  photon  density  rate  equations  are  shown  equivalent. 
The  mixed  flux-population  equations  are  valid  within  the  confines  of  the 
point  model. 


